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ABSTRACT'''' 

A  target  (deer)  is  moving  along --'a  path  in  the  forest. 

The  path  is  partly  obscured  from  an  observer  (hunter) 
by  trees  which  are  randomly .dispersed.  In  order  to 
detect  the  target  there  should  be  a  visible  window 
along  the  pat)j  of  length  ^lror  larger.  If  a  target  is 
detected,  the  hunter  starts  shooting  at  it. 
shooting  tri:al  requires  a  visibility  window 
of  length  /  .  The  hunter  continues  with  the 
trials  until  the  target  is  hit  or  disappears 
shadowed  portion  of  the  path.  The  present  paper  pro¬ 
vides  a  methodology  for  approximating  the  probability 
of  detecting  a  target,  and  lower  and  upper  bounds  for 
the  probability  distribution  of  the  total  number  of 
shooting  trials  along  the  visible  portions  of  the  path. 
Lower  and  upper  bounds  for  the  probability  of  hitting 
the  target  are  provided  too.  These  bounds  and  approxi¬ 
mations  are  derived  under  the  assumption  that  the 
obscuring  elements  constitute  a  Poisson  random  field.  C" 
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0.  Introduction 


A  hunter  is  trying  to  detect  and  hit  a  deer  in  a  forest. 

Suppose  that  a  deer  is  moving  along  a  path  in  the  forest  and 
the  hunter  is  located  among  the  trees  at  some  distance  from  the 
path.  The  path  is  only  partially  visible  to  the  hunter;  the 
invisible  (shadowed)  portion  of  the  path  is  obscured  by  the 
trees  which  are  dispersed  randomly  between  the  hunter  and  the 
path.  A  deer  can  be  detected  by  the  hunter  if  at  least  a  certain 
part  of  it  is  visible.  This  occurs  only  if  at  least  one  of  the 
visible  segments  of  the  path  is  sufficiently  long. 

After  detection  of  a  deer,  in  a  visible  segment,  the  hunter 
starts  shooting.  The  deer  continues,  however,  to  move  along  the 
path  in  the  same  pace.  During  each  shooting  trial  the  deer  crosses 
a  length  of  t  of  the  path.  Thus  the  number  of  shooting  trials  in 
each  visible  segment  depends  on  the  length  of  the  segment.  The 
shooting  trials  stop  when  the  deer  is  either  hit  or  enters  an 
invisible  portion  of  the  path.  When  the  deer  enters  another  visible 
segment,  it  has  to  be  detected  agin.  For  simplicity  we  assume 
that  the  shooting  trials  are  Bernoulli,  with  probability  of  failure 
q,  0<q<l . 

The  problem  of  deer  hunting  can  be  treated  as  a  two  dimensional 
shadowing  problem.  The  hunter  is  located  at  a  point  0  in  the 
plane,  the  deer  moves  along  a  curve  C  in  the  plane,  and  the  trunks 
of  trees  can  be  described  as  random  disks  dispersed  between  0  and  C. 
Two  dimensional  random  shadowing  problems  were  previously  studied 


by  Chernoff  and  Daly  [1].  Likhterov  and  Gurin  [2],  Yadin  and 
Zakcs  [3,4].  The  methodology  developed  in  the  present  paper  is 
also  applicable  to  three  dimensional  versions  of  the  above  problem. 
For  example,  if  a  hunter  tries  to  shoot  down  a  bird  whose  flying 
course  is  partially  obscured  by  crowns  of  trees.  The  three  dimen¬ 
sional  shadowing  problem  was  previously  studied  by  Yadin  and 
Zacks  [5] . 

In  the  present  study  we  develop  approximations  for  (a)  the 
probability  of  detection;  (b)  the  probability  distribution  of  the 
maximal  number  of  shooting  trials  N;  and  (c)  the  probability  of 
survival  of  the  deer  (bird) .  We  also  provide  numerical  examples 
to  illustrate  the  goodness  of  these  approximations. 

In  section  1  we  review  the  notions  and  definitions  concerning 
random  Poisson  shadowing  fields  in  the  plane,  measures  of  random 
visibility  on  star  shaped  curves,  and  introduce  the  notion  of 
x-reduced  measures  of  visibility.  This  notion  enables  us  to  define 
detection  probabilities  in  terms  of  the  distribution  of  the 
T-reduced  measure  of  visibility.  In  addition  we  can  also  arrive 
at  lower  bounds  to  the  number  of  Bernoulli  trials,  N,  along  the 
visible  portion  of  a  curve.  The  moment  generating  function  (MGF) 
of  the  x-reduced  measure  of  visibility  yields  lower  and  upper 
bounds  to  the  probabilities  of  a  successful  hunt.  Section  2  dis¬ 
cusses  the  moments  of  the  T-reduced  visibility  measure  and  intro¬ 
duce  an  approximation  to  its  distribution.  With  the  aid  of  this 
approximation  we  obtain  bounds  for  the  MGF  mentioned  above. 


Section  3  and  4  present  approximations  to  the  cumulative  dis¬ 
tribution  functions  of  the  T-reduced  visibility  measure  V(x),  and 
of  the  number  of  trials,  N,  along  the  visible  portions.  Section 
5  is  devoted  to  the  problem  of  determining  the  probability  of 
simultaneous  visibility  of  n  points  along  the  curve  C,  and  the 
recursive  determination  of  the  moments  of  V(t).  Section  6  provides 
some  special  cases  of  two-dimensional  models.  Finally,  in  Section 
7  we  provide  a  numerical  example.  The  methodology  presented  here 
for  two-dimensional  shadowing  problems  can  be  generalized  in  a 
straight  forward  manner  to  three-dimensional,  by  applying  similar 
modifications  to  the  results  of  [5]  . 

1 .  The  Model,  Measures  of  Visibility  and  Failure  Probabilities 
Suppose  that  the  hunter  is  located  at  the  origin,  0,  and  let 
C  denote  the  path  of  the  deer.  C  is  assumed  to  be  a  smooth  star 
shaped  curve,  defined  by  a  piece-wise  differentiable  function 
r(s),  sT  <  s  <  s  .  representing  the  distance  from  0  to  C  in 

±j  u  — 

orientation  s.  The  polar  coordinates  of  a  point  Pg  on  C  are 

(r(s),s).  The  end-points  of  C  are  P  and  P  .  The  length  of  C  is 

~SL  ~su 

s 

r  u 

(1.1)  L  =  I  £(s)ds 

SL 

where 

Us)  =  [r2  (s)  +  (  5-  r  ( s)  )2]1/2 


The  trees  in  the  forest  are  presented  by  random  disks  dispersed 
in  a  region  between  0  and  C.  Each  random  disk  is  characterized 


by  coordinates  (p,0,y),  where  (p,0)  are  the  polar  coordinates  of 
its  center  and  y  is  its  diameter.  The  coordinates  (p,S,y)  belong 
to  a  set  S  in  satisfying  conditions  which  assure  that  0  is  not 
covered  and  C  is  not  intersected  by  random  disks.  Let  8  be  the 
Borel  a-field  on  the  sample  space  S,  and  let  N{b)  designate  the 
number  of  disks  whose  coordinates  belong  to  a  set  B  of  B.  We 
assume  that,  for  each  Be8,  N{B}  is  a  random  variable  having  a 
Poisson  distribution  with  mean 


(1.2)  v{B>  =  A  jjj  H  (dp  ,d9  )  dG  (y  j  p ,  0  )  , 


where  G(y|p,0)  is  the  conditional  CDF  of  y,  given  (p,0),  and 
H(dp,d0)  is  a  a-finite  measure  of  (p,G).  Such  a  random  field 
of  disks  is  called  a  Poisson  random  field. 

A  point  Pg  on  C  is  said  to  be  visible  if  the  line  segment 

0Ps  is  not  intersected  by  any  random  disk.  A  point  which  is  not 
visible  is  in  a  shadow.  The  measure  of  total  visibility  on  C 
is  defined  as 


(1.3) 


C  U 

=  J  I(s)Ms)< 


where  I(s)=l  if  Pg  is  visible,  and  l(s)=0  otherwise.  Notice  that 
V  is  a  random  variable  representing  the  total  length  of  the 
visible  portion  of  C.  V  is  a  sum  of  a  random  number,  M,  of 
visible  segmentsof  C  having  random  length  X,,  X~,  ...,  XM;  i.e. 


A  target  is  detected  only  if  there  exists  at  least  one  visible 
segment  of  length  greater  than  the  minimal  path  length  xq  required 
for  identifying  the  target.  In  order  to  develop  a  formula  for 
the  probability  of  detecting  a  target,  we  introduce  the  notion 
of  r-reduced  visibility  measure,  V(t),  which  is  the  total  length 
of  visible  segments,  each  one  reduced  by  x  units,  i.e., 

M 

(1.5)  V ( x )  =  I  (X, -x) 

i=l  1 

where  a+  =  max(a,0).  The  probability  that  a  target  is  not  detected 
is 

(1.6)  P0(T0)  =  Pr{V(xQ)  =  0}. 


On  the  other  hand,  the  probability  that  C  is  completely  visible  is 


(1.7)  p.^  =  Pr(V(x)  *  L-x),  for  all  0<t<L. 


Indeed,  when  C  is  completely  visible,  M=1  and  X1=L.  Let  N  denote 
the  number  of  shooting  trials,  after  detecting  a  target.  If  a 
single  shooting  trial  requires  a  segment  of  length  x  to  be 
completely  visible,  then 


(1-8) 


N 


l<xi.-ToVTi 


/ 


where  [a]  is  the  maximal  integer  not  exceeding  a.  Notice  that 


K'.' 


h  ^  -  » 

P 


>- 


M 


M 


(1.9)  7  E  (X.-T  -x)  £  N  <  7  E  (X.-T) 

T  i=l  1  T  i=l  1 


Hence,  according  to  (1.5)  and  (1.9), 


(1.10) 


V(t1) /t 


1  N  1  V(T  )/x 


where  x .  =  t  +  t  , 
1  o 


If  the  probability  of  failure  in  each  shooting  trial  is  q,  and 
the  shooting  trials  are  independent  (Bernoulli) ,  the  number  of 
shooting  trials  required  until  the  first  success,  J,  is  distributed 


geometrically.  Accordingly,  the  probability  of  failure  (not 


N 


hitting  the  deer)  is  Q  =  E{q  }.  Thus,  according  to  (1.10),  lower 


and  upper  bounds  for  Q  are,  respectively,  Qq  and  Q.^  where 


(1.11) 


0i  =  E{qV(V/T) 


,  i  -  0,1 


Notice  that  is  the  value  of  the  MGF  of  V(x.)  at  the  point 
t  =  (log  q) /x  . 


2 .  The  Moments  and  Moment  Generating  Function  of  V(x). 

For  the  sake  of  determing  the  moments  of  V(x)  we  introduce 
the  following  definition  of  this  measure, 

s. 


(2.1) 


V(x) 


■  / 


U,T 


I  T(s)  £.(s)ds 


L,x 


JL  AJl i  jJL  nfn  It&tf  ‘iffrrT  if  m  Itlnr  jCwUtL Jfc JL  *•  .  I.  ,,*«> 


where  I ^.  ( s )  =  1  if  a  segment  of  C  of  length  t,  centered  at 

(r(s),s)  is  completely  visible,  and  I  (s)  =  0  otherwise.  sT 

and  sn  are  direction  coordinates  of  points  within  C,  of  distance 
t/2  along  C  from  sL  and  sy  respectively.  More  formally,  let 


(2.2)  L (s)  =  J  £(y)dy. 

SL 

Then,  s  =  L-1  (t/2)  and  sIT  =  L_1(L-x/2). 

J-j  /  i  U  /  T 

The  n-th  moments  of  V(t)  is  thus 

SU,T 

nn(x)  =  E{ (  j  IT (s) £{s)ds)n} 

SL,t 

(2.3) 


A 


n,  x 


n  n 

E{  TT  I  (S.  )  }  TT  l  (S.  )dS.  . 

i=l  T  1  i=l  ii 


The  set  A  is  the  simplex 
n ,  x  * 


(2.4) 


A 


n,  t 


s  );  sT  <s.< 
n  l,t—  l 


.  .  <s  <sTT  ; 
n—  U,  T 


n 

Furthermore,  E{  ir  I  (s.)}  is  the  probability  that  the  union  of 
i=l  1  1 

n  segments  of  C,  each  one  of  length  x,  centered  at  n  points 
having  direction  coordinates  s^<...<sn,  is  completely  visible. 

This  probability  is  designated  by  Pn(s]_'  •••>  sn;x).  Thus  the 


n-th  moment  of  V(x)  is 


k.' 


(2.5) 


Hn(t)  =  n!  J 


J 


Pn(sl 


n,  t 


TT  *-  (S  .  )  dS  .  . 
1=1  1  X 


The  method  for  determining  p  (s.,  s  ;x)  and  y  (t)  is  based 

n  in  n 

on  a  general  methodology  developed  by  Yadin  and  Zacks  [3,4]  for 
the  special  case  of  x=0.  The  main  results  of  those  papers,  with 
the  modifications  required  for  x>0,  are  presented  in  Sections 
4  and  5 . 


3 .  An  Approximation  to  the  CDF  of  V(x) 

The  cumulative  distribution  function  (CDF)  of  V(x)  is  a 
mixture  of  a  two-point  distribution  concentrated  on  {0,L-x}  and 
a  distribution  concentrated  on  the  interval  (0,L-x).  For  the 
purpose  of  presenting  the  approximation  discussed  below,  we 
consider  a  normalized  measure  of  visibility  W(x)  =  V(x)/(L-x), 
which  is  concentrated  on  [0,1].  The  CDF  of  W(x)  can  be  represented 
as 


!  0 


,  if  w<0 


(3.1)  Fx(w)  =  )  P0(x)  +  (l-p0(x)-Pl)  F*  ( w)  ,  0<w<l 


\  .i. 


,  1<  w 


If,  for  example,  G(y|p,0)  is  absolutely  continuous  then  F*  ( w)  is 
an  absolutely  continuous  CDF  on  (0,1).  Let'u  (x)  denote  the  n-th 
moment  of  W(x).  Obviously,  nn(t)  =  (L-x)nun(x)  ,  n=l,2,... 


Furthermore , 


I: 


(3.2)  un(T)  =  px  +  (l-po(r)-Pl)  J  w ndF* ( w) 

0 


n  =  1 ,  2 ,  .  .  . 


Applying  the  Dominated  Convergence  Theorem  one  immediately  proves 

that  lim  u  (x)  =  p.  for  all  x>0. 

n  1  — 

n-**» 


Explicit  expressions  for  po(x)  and  F*(w)  are  not  available. 

We  apply  here  a  beta  approximation  to  F*  (w )  and  provide  a  numerical 

approximation  to  pq(t).  This  type  of  mixed-beta  approximation  was 
applied  also  in  [3,4,5].  As  will  be  shown  in  Section  6,  in  some 
special  cases,  the  first  ten  moments  of  W(x)  and  of  the  mixed-beta 
approximation  are  very  close.  This  indicates  that  in  those  cases 
one  has  a  highly  effective  approximation.  In  cases  where  the  moments 
are  not  in  agreement  better  approximation  should  be  attempted. 

The  approximating  beta-mixture  CDF  is  given  by  the  formula 


(3.3) 


Ft(w) 


0 

<  Pq(t)  +  (l-P0(x)-Pl)Iw(aT,eT), 
1 


,  if  w<0 
0<w<l 

if  1<  w  , 


where  Iw(a,B)  ,  0  <_  w  <_  1 ,  0<a,8<°°,  denotes  the  incomplete  beta 
function  ratio.  The  probability  of  complete  visibility  of  C,  p^, 

is  determined  by  the  shadowing  model,  as  shown  later.  The  values 
of  P0(t) ,  aT  and  ST  are  determined  by  equating  the  formulae  of  the 

first  three  moments  of  F^ (w)  to  those  of  W(x).  As  shown  in  [3] 
these  approximating  parameters  are  given  by 


(3-4) 


a  = 

X 


2(y*(x))  -y*(x)  (y*(x)  +  y*(x)) 


y*  (x) y*  (x)  -  (y*  (x)  ) 


(y£(x)-y*(x))  (y*(T)-y*(x)) 


y*  (x)  y*  (x)  -  (y  £  (x) ) 


(3.5) 


where 


(3.6) 


PQ(x)  ■  l-p1-yj  (x)  (aT+BT) /oT 


Mj(T)  -  Un(T)-Pl 


n=l ,2,3  . 


Notice  that  y*(x)>0  for  all  n=l,2,...  and  that  y*(x)>y*  .(x). 

n  —  n  —  n+1 

As  revealed  in  formula  (3.4)  the  solutions  for  a  and  B  might 

x  x 

be  sensitive  to  the  numerical  accuracy  of  the  moments,  p  (x  )  is 

o  o 

the  mixed-beta  approximation  to  the  probability  of  no  detection, 

P  (x  )  . 
o  o 

Likhterov  and  Gurin  [2]  developed  a  theory  of  random  coverage 
on  the  real  line,  according  to  which  P0(tq)  can  be  obtained  by 
solving  certain  integral  equations.  Such  a  solution  requires 
generally  numerical  iterative  techniques.  The  method  of  approxi¬ 
mating  P0(t0)  from  the  moments  of  W(x)  requires  a  numerical 
approximation  too,  but  is  simpler  to  execute. 


4 .  Bounds  for  the  CDF  of  N  and  for  Q 

Inequality  (1.10)  yields  lower  and  upper  bounds  for  the  CDF 
of  N.  Indeed,  from  (1.10), 
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(4.1)  Ft  )<Pr{N<n}<FT  (y^-) 

Tq  L-Tq  -  ~  ~  T1 


The  CDF’s  in  (4.1)  can  be  approximated  by  the  mixed-beta  CDF 
(3.3).  According  to  (1.11),  the  lower  and  upper  bounds,  for 
the  failure  probability  Q,  are  the  value  of  the  MGF  of  W(t^)  i=0,l, 
at  the  point  t=^-(L-T^)  log  q.  Let  Gt  (t)  indicate  the  MGF  of  W(t)  . 

This  function  can  be  expressed  in  terms  of  the  moments  of  W(x)  as 


(4.2) 


.  °°  u*(x)  _ 

Gx(t)  =  1  +  pL(e  -1)  +  E  —  tn, 

n=l 


Since  y*(T)  4-  0  as  n  grows  the  infinite  series  in  (4.2)  converges 
faster  than  et,  and  therefore  a  small  number  of  terms  will  often 
provide  a  good  approximation.  Another  method  of  approximating 
Gx(t)  is  by  employing  the  MGF  of  the  mixed-beta  distribution  (3.3) 
with  p(x),  and  8^. 


5 .  Visibility  Probabilities  and  Moments  for  Two  Dimensional  Models 
In  Section  2  we  introduced  the  probabilities  Pn (si , • • • , sn ? t ) . 

In  the  present  section  we  discuss  the  methodology  for  determining 
these  functions  in  two  dimensional  models,  and  indicate  how  the 
moments  can  be  determined  from  (2.5)  in  a  recursive  fashion.  In 
addition  we  provide  an  explicit  formula  for  the  determination  of 
the  probability  of  complete  visibility,  p^.  We  start  with  the 
presentation  of  the  various  functions  required  for  the  determination 
of  Pn (s^ , • • • , sn;  0)  and  y  (0),  and  continue  with  the  modification 
needed  for  the  general  case. 
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Consider  the  star-shaped  curve  C,  specified  by 

(5.1)  C={ ( p  ,  0 ) ;  p=r ( 0 ) ,  sT<e<s  }  , 

where  r(0)  is  continuous  and  piece-wise  differentiable  on  [sT,s  ]. 

L  u 

We  assume  that  shadows  on  the  curve  C  are  cast  by  disks,  with 
centers  which  are  randomly  distributed  within  a  strip  S,  bounded 
by  the  curves 

(5.2)  U  «  {  ( p  ,  0 ) ;  p-u ( 0 ) ,  0  <9<0  } 

Lt —  U 

and 

(5.3)  W  =  {  ( p  ,  @ )  ;  p«w(0),  0. <0<6  } 

Li —  —  U 

where  ®l<sl<su<®u  an<^  u(0)<w(0).  For  simplicity  of  the  exposition 

we  consider  a  standard  Poisson  shadowing  process,  according  to 
which  the  centers  of  random  disks  are  uniformly  distributed  in  S 
and  their  diameters  are  distributed  independently  of  their  location, 
i.e. 

(5.4)  H (dp ,d@ ) dG (v | p , 0 )  —  pdpdedG(y). 

In  addition,  u(0),  w(0)  and  G(y)  should  satisfy  conditions 
which  insure  that  the  origin  0  is  uncovered  and  C  is  not  intersected 
by  random  disks.  In  particular,  these  conditions  require  that 


G(y)  will  be  concentrated  on  an  interval  [a,b]  and  that  the 
distance  of  U  from  0  and  of  W  from  C  would  be  at  least  b/2 


sa 


¥3W 


I 


(see  [4]).  For  the  purpose  of  obtaining  the  visibility  prob¬ 
abilities,  we  introduce  the  auxiliary  functions 

s  w(9 ) 

t  /  /  G  (y  (p , d-9 ) ) pdpd9  ,  t>0 

I  s-t  u(9) 


(5.5) 


K  (s,t)  = 


,  t<  0 


(5.6) 


K+(s,t)  = 


s+t  w(9) 

/  /  G (y (p , 9-s) ) pdodfi  ,  t>0 


s  u  ( 9 ) 


,  t<  0 


where 


(5.7) 


y (p , 9-s)  =  y (p , s-9)  = 


2psin | 9-s 


,  if  |  0  — s  |  <tt  / 2 


,  if  !  9  —  s  |  >tt  / 2 


Notice  that  y(p,9-s)  is  the  maximal  diameter  of  a  disk  centered 
at  (p,9)  which  does  not  intersect  the  ray  with  direction  s.  Thus, 
\K_(s,t)  and  XK+(s,t)  are  the  expected  number  of  disks  in  the  region 
S,  having  centers  with  direction  coordinates  in  [s-t,s]  and  [s,s+t], 
respectively,  which  do  not  intersect  the  line  segment  0Pg. 

Let  v{S}  denote  the  expected  number  of  disks  in  the  region 
S,  i.e. , 
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Wv:v.v •  ./.■..Vs-'v'* ■: 


(5.8) 


v 


{  S}=A  J 


J  pdpde  =jl  (w  (6)-u  (6))d0. 

U  (0  )  0 

Li 


Accordingly,  the  expected  number  of  disks  which  can  cast  shadows 
on  C  is 


(5.9) 


vx  =  v { S }-A [K_ 


|SL'SL  -6LI+K+|SU' 


6  -s  )  ] 
u  u 


•Thus  the  probability  of  complete  visibility  is  p^=e-vl.  The 

probability  p  (s. , . . . , s  ; 0)  fo  simultaneous  visibility  of  n  points, 

with  sT<s,<...<s  <s  is 
L—  1  n—  u 


(5.10)  Pn(si, . . . ,sn;0)  =  exp{-v{S}}exp{AK_(s1,s1-0L) 


n-1 


+AK 


+  <VW  +  X  2 


si+1  si)+K_(3.+1,ii±i-ii)]} 


Finally,  we  present  the  recursive  equations  for  the  calculation 
of  the  moments  y  (0)  .  Define  the  sequence  of  functions 

(5.11)  <i>Q(s)  =  exP  ( AK_ (s, s-sL)  } 

and  for  each  j>^l 

s 

~  l  f  +  xK+(y,^-ndy. 


(5.12) 


Then,  for  each  n>l,  the  moments  of  W(0)  are 


(5.13) 


un(°) 


ru 

=  jr  exp{-v{S}}  J  £(s)^n_1 (s)exp{ AK+(s, 0u-s) }ds, 


In  order  to  generalize  the  above  results,  one  has  to  adjust  the 
functions  K_(s,t)  and  K+(s,t)  in  the  following  manner.  For  a 
given  orientation  s,  define 


(5.15) 


t_  ( t , s )  =  s-L  iL(s)-t/2) 


(5.16)  t+(T,s)  *  L_1(L(s)+t/2)-s. 


The  function  p  (s.  ,...,s  ;x)  is  obtained  bv  substituting  in  (5.12) 
n  i  n 

K_(s-t_(x,s) ,t-t_(x,s) )  and  K+ (s+t+ (x , s) ,  t-t+(x,s))  for  K_(s,t) 
and  K+(s,t),  respectively.  Similarly,  the  moments  un(x)  are 
obtained  by  replacing  K_(s,t)  and  K+(s,t)  in  formulae  (5.12)  - 
(5.13)  by  K_(s-t_(x,s) ,t-t_(x,s) )  and  K+ (s+t+ ( x ,s) , t-t+ (x , s) ) , 
respectively;  and  replacing  L  by  L-x. 


.  Some  Soecial  Cases  of  Two  Dimensional  Models 


6.1  Circular  Path  and  Annular  Region 


The  curve  C  and  the  region  S  are  specified  by  the  functions 
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u (0)  =  u, 


-■s-<eT<8<0  <  tt/2 

Z  —  Li  —  U  — 


w  (9 )  =  w 


-^<eT<0<6  <  tt/2 

2  —  L  —  —  n  — 


r  (s)  =  r. 


<  sL  £  »  i  su  <  eu  . 


Furthermore,  we  require  that 

0  <  b/2  <  u  <  w  <  r-b/2 

Due  to  the  symmetry  on  the  circle,  one  can  show  (see  [4])  that 
the  K-f unctions  in  the  annular  case  do  not  depend  on  s,  and  satisfy 


(6.1)  K+(s,t)  =  K_ (s , t)  *  K*(t) 


Explicit  expressions  for  K*(t),  in  the  case  of  a  uniform  distri¬ 
bution  G(y)  on  [a,b],  was  derived  in  [4],  Let  K*(t)  =  K  (t)-K  (t) 

w  u 

The  functions  K  (t)  and  K  (t)  can  be  obtained  from 
w  u 

/  0  ,  if  t  <  sin-1(-^) 


Kv(t)  = 


(t) 


K<2>  (t) 


,  if  sin'1  (^  £  t  <  sin"1(^) 


,  if  sin'1^)  1  t  <  tt/2 


K(3) (t) 
v 


,  if  t  >  tt/2  , 


where 


-2v  cos ( t)  ] 


k'11  ,t)  =5^  if  [Hv2^2)1'2 


a  2,.  -1  .  a.,.  a  r  2  2.1/2 

-v  [ t-sin  <2v)l+  y^tU'v  _a  ) 


a  cotan  ( t)  ] } , 


K^2)(t)  =  K^1)  (sin"1(^))+  (t-sin-1(|^)) 

2  2 

a  +  ab+b  ,1,,  2,2,  1/2  t  , 

- -  [_(4V  ~b  '  “  cotan  ( t )  ]  , 


and 


k<3»  (t)  -  I-421  (f)  *(t-  |)  -  afiablfe!) 


Due  to  the  simple  and  symmetric  nature  of  C  and  S,  the  adjustment 
factors  t_(t,s)  and  t  (t,s)  satisfy 


(6.2)  t_ (i|S)  =  t+(T,s)  =  t/2r 
Accordingly,  formulae  (5. 8) -(5. 10)  are  replaced  by 

(6.3)  v{S}  *  ^(w2-u2) (0u-0L) 

(6.4)  \>1  -  v{S}  -  X[K*(sL-9L)+K*(6u-su)]  , 


and 

pn*sl' *  * ' 'sn?T*  =  exp{-v{S}}exp{X[K*(s1-6L-  |^) 

(6.5) 

T  n_1  Si+l“Si  T 

+  K*<9u-Sn  -  h)]  +2X  K*(^r-^  -  h  »• 


-17- 


■l*  ■  i  a  •_*  tLa  " 


For,  determining  the  moments  u  (x)  we  compute  recursively 


(6.6)  ^Q(s)  =  exp{ AK* (s-sL - )}  , 


and  for  j^l  , 


(6.7) 


^ j  (s) 


exp{2AK*(^ 


dy' 


Using  these  functions  we  obtain , 


(6.8) 


pn(t)  - 


nrexpl- 


L-T~ ~} }'  1  Vl(y)eXpUK*(Vy-  27)}dy 


6 . 2  Linear  Path  C  and  Trapezoidal  Region  S. 

We  consider  here  the  case  where  both  U,  W  and  C  are  parallel 
lines  of  distance  u,  w  and  r  from  0.  Furthermore, 


(6.9) 


_  b  b 

0<-j£u<w<r--2 


In  this  case  it  is  simple  to  express  the  formulae  in  terms  of 
Cartesian  coordinates.  Accordingly,  a  point  P  on  C  has  the 
coordinates  (s,r).  Thus 
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«.  'V*'. 


(6.10) 


C  =  Us,r)  ;  sL  <  s  <  su)  , 

(6.11)  u  =  (e,u)  ;  eL  <  e  <  eu 

and 

(6.12)  W  =  {  (0  ,u)  ;  0.  <  0  <  0,  >  • 

It  was  shown  in  [5]  that  K_ (s , t) *K+ (s , t) =K (s, t) .  A  formula  for 
K(S/t)  is  given  in  [5].  We  adopt  the  convention  that  K(s,t)=0 
for  all  t<0.  For  the  t-reduced  case  we  obtain 

(6.13)  t_(T,S)=t+(T,S)=T/2 

Formulae  (5. 8) -(5. 10)  assume  in  the  t-reduced  case  the  form: 

(6.14)  v(S}  =  ^  (0u-0L) (w2-u2) 

(6.15)  vx  =  v{S}  -  X [K(sl,sl-0l)  +  K(sn'eu~su) 1  ' 

and 

(6.16)  pn(sL . .  *  exp{-v{S}} 

exp{X [K ( s1 , s1— 0L-  L)  +K(sn,0n-sn  -  ^)]+ 

X  l  [K(s.,s.-i)  +  K (s .  . , s .  -  ^)]} 


c 


where 


(6.17)  ii=rtan  (■j(tan~1  (^)  +  tan  1  (  ■i*1)  )  ) 

For  further  explanation,  see  [5] .  Finally,  determine  recursively 
the  functions 


(6.18)  ^o(s)  =  ex?{ ^K(s,s-6^-  j)  } 


and  for  j >1 , 


(6.19)  ^j(s)  =  f  4<j_1(y)exp{2\K(y,s(y)-  J) }  dy, 


where 


(6.20)  s  (y)  =rtan  (-^(tan  1  (^)  +  tan  1  (^) )  ) 


The  moments  of  W(t)  are  given  then  by 


s 

ru 

(6.21)  Mn(t)=  j^-exp{-v{S}  }  J  'i'n_1  (y)exp{  AK(y,  0u-y-  j) 


dy. 


7.2  Numerical  Example 

In  the  present  section  we  provide  an  example  which  demonstrates 
numerically  the  results  of  the  present  paper.  We  consider  the 
case  of  an  arc  C  and  annular  strip  S,  which  was  discussed  in 
Section  6.1.  The  parameters  of  this  case  are: 

0.  *  -tt/2,  st  »  -it/3,  s„  =  tt/3 ,  9  =  tt/2,  r*l,  w=. 6 ,  u=.  4 ,  X=5 . 

Lt  i-i  U  u 
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In  addition,  the  diameters  are  uniformly  distributed  over  the 
interval  ( .  1 ,  . 5 ) . 


In  Table  7.1  we  present  the  first  10  moments  of  W(t) ,  for 
x=0(.l).4,  which  have  been  determined  according  to  (6.6)— (6.8) - 
The  corresponding  moments  of  the  mixed-beta  distribution  (3.3) 
are  also  given  for  comparison. 


Vn 

t\. 

1 

2 

3 

4 

5 

6 

--I 

7 

8 

9 

10 

0.0 

.738 

.738 

.600 

.600 

.517 

.517 

.462 

.463 

.425 

.425 

.398 

.398 

.378 

.378 

.363 

.363 

.351 

.351 

■HE 

m 

0.1 

.704 

.704 

.561 

.561 

.479 

.479 

.427 

.427 

.393 

.393 

.369 

.369 

m 

.338 

.338 

.329 

.329 

m 

0.2 

Ml 

.526 

.526 

.447 

.447 

.399 

.399 

.368 

.369 

.347 

.347 

.332 

.332 

m 

.313 

.313 

.307 

.307 

0.3 

.641 

.641 

.497 

.497 

.421 

.421 

B9 

.349 

.349 

m 

.318 

.318 

am 

Ultrii 

BBitl 

mm 

MR 

Bwih 

■mm 

BKliM 

mm 

1  j 

1111™ 

0.4 

.614 

.614 

.471 

.471 

.399 

.399 

.359 

.359 

.334 

.334 

m 

.299 

.300 

.294 

.294 

.289 

.290 

TABLE  7.1  Moments  of  W(x)  (upper  line)  and  of  Ft (w)  (lower 
line)  for  x  =  0  { . 1) . 4  and  n  =  1 ,  . . . ,  10 . 


As  shown  in  Table  7.1,  the  first  ten  moments  obtained  from  the 
mixed-beta  CDF,  Ft<w),  differ  from  those  obtained  by  the  recursive 

formulae  only  at  the  4th  decimal  place.  This  reveals  an  excellent 
approximation  to  the  CDF  of  W(t)  by  F (w) ,  in  the  case  under 
consideration.  In  Table  7.2  we  provide  the  parameters  of  the 
mixed-beta  distributions  associated  with  Table  7.1. 


FY  srv  <r.  -  .  *r.  W.  w-w  «r.  «r_  w.  v.  .  .  v.  -  .  -  .  -  .  *.-  .  -  .  -  .  -  .  -  .  -  .  •'.  -  .  -  .  • 


T 

a 

P0<t> 

0 

.2353 

.0064 

.27 

3.3905 

1.8888 

.1 

.2559 

.0119 

.27 

3.0675 

2.0334 

.2 

.2747 

.0194 

.27 

2.8000 

2.1640 

.3 

.2917 

.0298 

.27 

2.6076 

2.3093 

.4 

.3069 

.0431 

.27 

2.4814 

2.4808 

TABLE  7.2.  The  Parameters  of  the  Mixed-Beta  Distribution 

F  (w)  for  t  =  0(.1).4.  (a  denotes  the  standard 
deviations . ) 


The  values  of  pq(t)  in  Table  7.2,  provide  the  mixed-beta 
approximations  to  the  probabilities  PQ(^0)  of  not  detecting  a 
a  target.  This  is  obviously  an  increasing  function  of  tq. 

Thus,  in  the  present  example,  if  tq  =  .1,  P0(T0)  =  .012  while  if 
tq  =  .4,  P0(T0)  =  .043.  p^  =  .27  is  the  probability  of  complete 
visibility  along  the  path.  Since  the  moments  of  the  mixed-beta 
distributions  (w)  fitted  so  well  those  of  W(t),  we  replace 

F  (T  Tt?  )  with  F_  (=-T2  -)  ,  i  =  0,  1.  In  Table  7 . 3  we  present 
Ti  L"Ti  Ti  L~Ti 

Fx  for  T.  =  0  ( . 1 ) . 4 ,  t  =  .1. 

Ti  L'Ti  1 


»n  aL  1  i 

The  values  of  Q^=E{exp{ t^W(T^) } }  where  t^  =  - - -  log(q) 

with  q=.8,  are  also  given  in  Table  7.3. 
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n^ 

— 

0.0 

0.1 

0.2 

0.3 

0.4 

0 

0.0064 

0.0119 

0.0194 

0.0298 

0.0431 

1 

0.0065 

0.0122 

0.0203 

0.0318 

0.0466 

2 

0.0075 

0.0147 

0.0255 

0.0411 

0.0613 

3 

0.0104 

0.0211 

0.0375 

0.0804 

0.0894 

4 

0.0166 

0.0332 

0.0577 

0.0906 

0.1307 

5 

0.0273 

0.0520 

0.0870 

0.1314 

0.1836 

6 

0.0435 

0.0785 

0.1253 

0.1819 

0.2458 

7 

0.0662 

0.1131 

0.1723 

0.2405 

0.3145 

8 

0.0960 

0.1557 

0.2269 

0.3052 

0.3865 

9 

0.1333 

0.2059 

0.2879 

0.3736 

0.4585 

10 

0.1782 

0.2628 

0.3532 

0.4430 

0.5272 

11 

0.2303 

0.3252 

0.4209 

0.5106 

0.5894 

12 

0.2890 

0.3914 

0.4884 

0.5735 

0.6423 

13 

0.3531 

0.4592 

0.5529 

0.6288 

0.6836 

14 

0.4208 

0.5261 

0.6115 

0.6738 

0.7117 

15 

0.4902 

0.5891 

0.6613 

0.7063 

0.7265 

16 

0.5582 

0.6448 

0.6992 

0.7249 

1.0000 

17 

0.6213 

0.6896 

0.7226 

1.0000 

1.0000 

18 

0.8750 

0.7183 

1.0000 

1.0000 

1.0000 

19 

0.7139 

1.0000 

1.0000 

1.0000 

1.0000 

20 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

Qi 

0.0704 

0.0967 

0.1273 

0.1621 

0.2000 

TABLE  7.3.  The  CDF  (7-^— -)  ,  with  T  ■  .1,  T.=0(.l).4, 

-  T.  L-T.  -  1 


L=su~sl;  and  the  corresponding  MGF  Q^. 


As  we  see  in  Table  7.3,  if  t=.1  and  t  =  .1  the  lower  bound 

o 

of  Q  is  .0967  and  the  upper  bound  for  Q  is  .1273.  If  however, 
To=0  then  .0704  <  Q  <  .0967. 


The  bounds  for  the  CDF  of  N  are  read  from  Table  7.3  in  a 

similar  manner.  For  example,  if  tq=0,  t^=.1+t  =.l  then  for 

n=6,  .0435  <  P{N  <  6}  <  .0785.  If,  T  =  .1  then  x,  «  .1  +  x  = 
—  —  —  o  i  o 

and  .0785  £  P{N  £  6}  £  .1253.  Thus,  from  the  first  two  columns 

of  Table  7.3  we  obtain  that,  when  tq  =  0,  the  expected  number 

of  trials,  E{N},  is  between  13.7  and  15.1. 
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